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^ ! Abstract 

(NJ ! It can be shown using operator techniques that the non- Hermitian PT- 

symmetric quantum mechanical Hamiltonian with a "wrong-sign" quartic 
\q ' potential — gx 4 is equivalent to a Hermitian Hamiltonian with a positive 

quartic potential together with a linear term. A naive derivation of the same 
i-g ■ result in the path-integral approach misses this linear term. In a recent pa- 

per by Bender et al. it was pointed out that this term was in the nature 
of a parity anomaly and a more careful, discretized treatment of the path 
integral appeared to reproduce it successfully. However, on re-examination 
of this derivation we find that a yet more careful treatment is necessary, 
^ ■ keeping terms that were ignored in that paper. An alternative, much sim- 

pler derivation is given using the additional potential that has been shown 
to appear whenever a change of variables to curvilinear coordinates is made 
in a functional integral. 

PACS numbers: 11.30.Er, 03.65.Db, 02.70.Dh 



1 Introduction 

Hamiltonians of the form 

H = ~p* - g(ix) N , (1) 



which are are PT-symmetric, but not Hermitian, have been shown ^ |2j to 
have a real, positive spectrum for N > 2. For N < 4 the energy eigenvalue 
problem can be posed on the real x-axis, but for N > 4 it must instead 
be imposed in the complex x-plane, along a contour which asymptotically 
remains within the Stokes wedges of the Hamiltonian. Such is indeed the 
case for N = 4, when the potential on the real axis is an "upside down" 
quartic. To that end we rewrite Eq. for iV = 4 in terms of the complex 
variable z: 

E = \vl~ 9Z 4 - (2) 

A general framework encompassing PT-symmetric Hamiltonians was given 
by Mostafazadeh t 3 , who showed that such Hamiltonians were related by a 
similarity transformation to an equivalent Hermitian Hamiltonian h possess- 
ing the same spectrum. In only a few cases is it possible to construct h 
explicitly, but this was done recently for Eq. (j2J) in Ref. [I], rediscovering 
an earlier result of Buslaev and GrecchijS]. The calculation was performed 
using the particular real parametrization 

z = -2iVl + ix (3) 

of the contour in the complex plane on which Eq. @ was defined, and using 
operator techniques to perform the similarity transformation to h, whose 
form is (a = 16g) 

We should emphasize that this latter form is obtained only after a Fourier 
transform, which means that y is really a momentum rather than a position, 
as will become apparent in the path- integral treatment. When a proper 
dimensional analysis is made|6j, the linear term turns out to be proportional 
to h, which means that it is in the nature of an anomaly. 

In fact we will present various path-integral treatments, of various levels of 
sophistication. In the next section we give the simplest of these, a continuum 
version, which fails to reproduce the anomaly. Realizing that for a general 
change of variables in a functional integral there is an additional term, AV, to 
be included over and above the functional determinant, we show, in Sectional 
how AV is traded for the anomaly in an intermediate functional integration. 
In Section 4 we give the discretized version of this calculation, in order to 
provide a benchmark for comparison with the results obtained in Section 
5 by two variants of the discretized calculation of Ref. 0. In effect this 
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latter section gives an alternative derivation of AV, as a consequence of 
the mismatch between the functional determinant and the coefficient of the 
kinetic term when our particular change of variable is made. 



2 Naive Continuum Treatment 



The original path integral is the single Euclidean functional integral 



Z = J [Lhflexpl- J dt 



(5) 



where the subscript C is to remind us that the integrals are to be performed 
on an appropriate curve in the complex ip plane. Making the change of 
variable 



7p = -2iy/(l+i(p), 
in analogy with Eq. (j3*|) . we obtain 



Z 



[Dp] 



exp 



dt 



1 2 



2 1 + i<p 



-a(l + i<pf 



(6) 



(7) 



where the <p integrals are now along the real axis. We now rewrite this as a 
double functional integral in tp and tt by means of the identity 



1 



DetV(l+*» 



1 

N 



[Dtt] exp 



dt-(l + i(p) ( tt - % — - — 



where N is an appropriate normalization constant. This leads to the expres- 
sion 



Z 



1 

N 



[Dip] [Dtt] exp 



dt 



+ i p)tt 2 — iTfp — ct(l + i ipY 



(9) 



Replacing mp by -Trip under the t integration, we have an exponent that is 
quadratic in <p, so the <p functional integration can be performed. The result, 
after applying the compensating rescalings <p> — > <p/\f2a and ix — ► n\/2a, is 



[Dn] exp 



dt 



1 



-7T — v2tt7r ( 1 — —TT ) + —IT 



a 



(10) 



The middle term in the integrand is a perfect derivative, and so may be 
discarded, leaving the equivalent Lagrangian I in the exponent, written in 
terms of tt: 



[Dtt} exp 



dt 



1 - 2 , " 4 

— TT H TT 



ID 



except that the linear term is missing. The simple classical calculation we 
have just performed is unable to obtain this term. Note that the Fourier 
transform of the Schrodinger treatment occurs here naturally, since Z is 
expressed in terms of the momentum variable tt. 



3 Correct Continuum Treatment 



It is shown in various standard books on functional integration, for example 
[7JIB], that when a general change of variables, such as that of Eq. (JBJ), is made, 
an additional potential term AV must be included in the Lagrangian. Ulti- 
mately this term, which is actually of order % , is derived from the discretized 
form of the functional integral when the particular form of discretization 



a 

<Pn = o(y?n+l + V?n) 



(12) 
(13) 



is adopted, where a is the lattice spacing, corresponding to Weyl ordering 
in the operator treatment. The importance of this prescription was indeed 
emphasized in Ref. [Hj. 

The general form of AV for a change of variable from ip to ip is [7] 



AV = - 



1 
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d(p 





which for the particular transformation of Eq. (jHJ) turns out to be 

AV=-i_U 
32 1 + iip 

The correct version of Eq. (JJJ) is thus 



(14) 



(15) 



[Dip] 



Dety/(1 +i<p) 



exp l—fdt 



1 1 



21 + iip 321 + ^ 



a(l + iipf 



It is now possible to write down a variant of the Gaussian identity 

1 x 2- 



1 



Dety/(1 +i<p) N 



[Dn]exp{- jdt\{l + iip) 



(16) 



(17) 



which serves to cancel the AV term as well as the kinetic term. In turn it 
introduces two additional terms in the exponent: (i) a term in i<p/(l + if), 
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which is a perfect derivative and so can be discarded under the t integration, 
and (ii) the anomaly J dt tt/4. 

Finally, after rescaling as before, the corrected version of Eq. (Illj) is 



[Dir] exp 



dt 



-71 

2 



. a 4 

— 7T H 7T 



4 Discretized Version 

In this section we will go through the discretized version of the previous 
calculation, in order to provide a standard discretized formula with which we 
can compare the results of the (corrected) calculation of Ref. jHJ and another 
calculation whereby the kinetic term is expanded in powers of the lattice 
spacing a. 

In place of Eq. (fTBj) we have 



n 



dip n 



== exp < —a 



1 vl 



1 



1 



V 1 + i( Pn 

The Gaussian identity we will use is 

1 



2 1 + iif n 32 1 + iip n 



- a{\ + ip n ) 2 



.(19) 



^ J d7r n e-^- B ^ 



VI + Wn N 

where A = a(l + i(p n ) and B = {i(p n + j)/(l + %(p r 
Written out in full this is 



(20) 



1 



V 1 + Wn 



^ J dn n exp 



1_ 2 . _ . _ 1_ 

-7T n (l + Up n ) - l7l n Lp n - ~TX n 



+ 



1 iif> r 



Pi 



+ 



41 + iip n 21 + iip n 32 1 + iip % 
Neglecting the term in i(p n j (1 + i<p n ) because the identity 



log 



1 + i(f n+ l 



i(up n 



+ 0(a 3 



(21) 



1 + i(p n J l + iip. 

shows it to be a perfect difference up to a correction of order a 3 , we obtain 
1 



Z = — n J dir n dip n exp 



-(1 + itPn)^ - ilT n (f n - ^7T n - a(l + i(f n ) 2 ) 



[22) 
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Thus the AV term has been cancelled, and we are left with the anomaly 
|a7r n in the exponent. Now 7r„0 n + (p n ^n is a perfect difference: 



(23) 



So now we can "integrate by parts", in the form n n (p n — > — (p n n n . 
Changing the integration measure from J d7t n d(p n to / d7T n d(p n , which does 
not introduce any additional factors, Eq. (J22J) becomes 



Z = — Y[ J dir n d<p n exp J -a a ( (ip n - i) + + 27r n )) 

Tl \ L 

1 /_ 2 \ 2 1 _ 1 1 

16a V n > 4 J J ' 

having dropped a perfect difference proportional to 7r n . 

Now we rescale: <£? n — > (p n /^/2a and 7r n — > n n ^/2a and perform the (p n 
integration, with the result 



z = n / rf7r « ex p 



a 



/a 



4< + \/ 2 7r "^ n 



2*" 



a _ 

— 7T r , 



This is the desired result, provided that we can neglect the term 7r^7r n . The 
identity 



3a7r^7r n 



7T 



n+l 



n 3 n - (7T n+ l - 7T n ) 3 /4 



(24) 



shows that is a perfect difference up to a correction of order a 3 , so that it 
can indeed be neglected. The resulting expression for Z is 

Z = ~[[J dir n exp j-a ^ - ^7r n + ^ j , (25) 

the discrete version of Eq. (|16p. Equation ()22|). from which this was derived, 
will be a reference point for the calculations of the following two sections, 
which arise out of the treatment of Ref . [0] . 

The basis of that treatment was an exact, discretized, treatment of the 
kinetic term. In the original functional integral written in terms of ip, the time 
derivative is defined as tp n = (tp n +i ~ ipn) / a - In terms of ijj n = —2iy/(l + iip n ) , 
this becomes ip n = tp n /A n , where 



A n 



1 + i<-p 



n+l 



ILfr. 



The relation between A n and y/1 + itp n can be written exactly as 



1 + 



ia(p r 



V 1 + i( Pn An \ 4A n \A + i( Pn 



(26) 



(27) 
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In this formulation the anomaly arises because the denominator A n in the 
expression for (fi n is not quite the same as that in the determinant of Eq. (}2*Tj) . 
That is, we have the expression exp[—a(p n 2 /(2A 2 l )]/V(l +iip n ), and we need 
to expand one of these in terms of the other. 



5 Expanding the Kinetic Term 



In this case we take A = a(l + i<p n ), x = n n and B = iip n /(l + i<p n ) in the 
Gaussian identity 



1 



, : — — f dxexpl — \(x — B) 2 \ 

1 f \ 1 . ,_ 2 1 

= J dir n exp <^ - -o(l + if n )^ n + %a ^Pn^n + -- 



(28) 



1 aft 



■2 



+ i(Pn 



Now we need to expand the kinetic term — ^aip n 2 /A\ in terms of \\B 2 
namely \\B 2 = —^a(p n 2 /(l + iip n ). First write Eq. (J27|) in terms of B: 

1 1 / aB I ; — 

— I 1 + -r-^V 1 + %LPn 



a/1 + Wn A. 



AAr. 



In terms of R = y/1 + i<p n / A n , this reads 

i = 1 + \aBR, 

which is a quadratic equation for R, with solution 

2 



(29) 



R 



l + VTTaB' 



So the kinetic term is -\atp n 2 /A 2 n = \\B 2 R 2 = \\B 2 + \\B 2 {R 2 - 1). 
Hence 



l-]a\B 3 + ^-(5\B i + \ 2 B e ) + ... 
4 32 



(30) 



Now we need the Gaussian identities, under Jdt exp{ — |A(x — B) 2 }: 

3;r 
T' 

Qx 2 3 

X + A 2 ' 
15x 4 A5x 2 15 

^T + ^~ A^' 



B 3 


= x 3 


B 4 


= x A 


B 6 


EE X 6 



(31) 
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to get 

e -l a ^/Al = e \XB* L + ^ + 0(q2) \ (32) 

Note that it is crucial to keep terms in Eq. (j3*Uj) which are nominally of order 
a 2 (and possibly higher). The point is that the Gaussian identities bring in 
terms of order 1/A, 1/A 2 etc., which means that such terms may actually be of 
order a. The net result so far is that we appear to have produced the anomaly, 
but with the wrong coefficient, 3/4 versus 1/4. However, it is important to 
realize that after these transformations we are left with — |a(l + i<p n )7t 2 in 
the exponent, as opposed to the — |a(l + i(p n )n 2 of Eq. (J22)- The difference 
between them is 

~a(l + ifn)^l = -\a(l + iy n )nl + Ua 2 <pnK 2 n 
= --a(l + i<p n )Trl + -XuBttI 

In order to implement Gaussian identities with the new kinetic term we need 
to write the correction in terms of a new A and B, namely A = a(l + i<p n ) 
and B = itp n / (1 + i(p n ), with the same property that XB = XB = iaip n . The 
exponential of the additional term ^XaBx 2 expands to 

exp^AcLBx 2 ) = 1 + -XaBx 2 + —X 2 B 2 x 4 + . . . = 1 - -ax + 0(a 2 ), 
4 4 32 2 

using the further Gaussian equivalences 

BX 2 EE X 3 -^, 

A 

Qr 4 12r 2 

BV = x e -^ + ^r. (33) 
A A 

Thus, the additional term precisely corrects the coefficient of the anomaly 
from 3/4 in Eq. (}3*2~|) to 1/4. A necessary ingredient for this to work is the 
lack of a term in 1/A 3 in the Gaussian equivalence for B 2 x A . 



6 Expanding the Determinant 

In this case we take x = n n , X = aA^ and B = itp n jA 2 n in the Gaussian 
identity 

X = ^/d*exp{-iA(*-fl)°} (34) 



1 f i f 1 ,2-2 • • a <Pl 1 

— J dx exp <--aA n n n + ia<p n n n + — | 



(35) 
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In fact we can solve Eq. (|2*7j) for a/1 + i(p n in terms of A r 



1 

Ar 



V 1 + Wn 



a/ 1 + Wn 



1 - -aB 
4 



Thus, in terms of the new B, R = ^1 + iip n /A n = 1 — aB /A. 
So the determinant can be written as 



= = — 1 + -aB + — a 2 B 2 + . 



16 



(36) 



(37) 



Now we use two further identities under the Gaussian integration of Eq. (J33J): 



B = x 
B 2 = x 2 - 



to obtain 



1 

A 



(39) 



The second term in ( ) correctly gives the anomaly when elevated to the 
exponent, but the third term, which was missed in Ref. jSj, appears to be an 
unwanted addition. The final result is 



Z = J| J dip n dtr n exp 



1 42-2 - • 1- 1 / 1 ■ 

jVn - Z7r "^" - 4^ + - a(l + If, 



(40) 

This differs from Eq. (|22|) in two ways: we have a term resembling AV, 
and the coefficient of 7f^ is rather than ^(l+i(p n ). The difference between 
and (p n in the interaction term is not important. 
Now 

aA 2 n 7? n = a(l + i^p n )^n + ^ 2 (l - R : 
So, since iaip n = ia{(p n — \<upn) = iafin ~ |AaS, 
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2 n n 2 

The additional terms expand to 



n 2 



-2 



4 nv ; 2 



i( a S) _ —(aB) 2 
2 V ; 16 V 7 
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and their exponential to 1 + (aB) 2 Xit 2 /32 + . . .. 

Again, because we have changed the coefficient of the kinetic term to 
|(1 + iifn) we strictly need to rewrite XB 2 in terms of A and B according to 
XB 2 = (XB 2 ) x (A/A), where the correction factor is 

A 1 + i(p n 1 i ^^ 5 s 2 



A Al 



l + 0(aBY 



The additional term does not in fact contribute to order a. Thus we need 
the final equivalence 

B 2 x 2 = * A --^ + Y 2 (41) 

for B 2 x 2 under Gaussian integration, which gives 1 + a 2 /(16A) + 0(a 3 ). 
When exponentiated this term precisely cancels the Al/-like term a/(16v4 2 ) 
in Eq. (gni) up to 0(a 2 ). 



7 Conclusions 

The naive functional integral formulation of the operator calculation can be 
regarded as the classical result, in that it takes no account of operator order- 
ing, or equivalently of any particular discretization, and so fails to produce 
the linear, anomalous term, which is of order H. We have shown an elegant 
method of producing this term in the continuous formalism, provided we 
take account of the additional AV^ term that has been shown (by careful 
discretization) to occur whenever a general change of variables is made. 

The calculation of Ref. jHj indeed attempted a careful discretization, but 
missed terms that were nominally of higher order in the lattice spacing a, but 
were in fact of the same order as the terms kept because of the particular na- 
ture of the Gaussian equivalences, which bring in factors of 0(1/ X) = 0(1/ a). 
Starting from the recognition that the coefficient of the kinetic term does not 
exactly match the argument of the functional determinant we proceeded in 
two ways, expanding either the kinetic term or the functional determinant. In 
both cases we used the discretized version of the corrected continuum calcu- 
lation as a canonical form with which to compare our results. The expansion 
of the kinetic term turned out to be the easier, needing only a change from 
the point variable (p n to the averaged variable <p n to obtain agreement. The 
method of expanding the determinant used in Ref. [0] eventually produced 
the same result after more lengthy manipulations, which incidentally showed 
that the simple passage from Eq. (57) to Eq. (58) in Ref. was not correct. 
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One can ask whether the same phenomenon of trading a. AV term for an 
anomaly by means of some analogue of Eq. (JT7j) is possible for other changes 
of variable. The answer seems to be in the negative: the parametrization of 
Eq. (JHJ) seems to be very special, as indeed it is in other respects. 

In general, the coefficient of the kinetic term in ip starting from a standard 
kinetic term in ifj is ~M, where 



to be compared with the expression for AV in Eq. (|14jl. For the two to be 
proportional we require 



or equivalently d 2 ip/dip 2 = const. This leads back to a relation of the general 
form of Eq. © up to a shift in ip. 

The main motivation for reformulating the quantum mechanical problem 
in path-integral terms in Ref. jH] was in order to attempt a generalization 
to higher dimensions, where there are indications from the non-Hermitian 
formulation that the theory is asymptotically free and naturally possesses 
a non-vanishing vacuum expectation value, making it an attractive alterna- 
tive to the standard Higgs model. It is therefore natural to ask whether 
additional potentials of the nature of AV, and the consequent production 
of an anomaly, occur in dimensions greater than 1. The answer appears 
to be negative, at least in the context of dimensional regularization. This 
is because algebraic field transformations lead to 0(h 2 ) terms proportional 
to (5^(0)) 2 in n spatial dimensions, which vanish in such a regularization 
scheme. This is consistent with a diagrammatic expansion, for which no ad- 
ditional terms are necessary jH]. The reason for their presence in quantum 
mechanics versus their absence in field theory has been analyzed in some 
detail by SalomonsonfTU]. 

Even though it appears that this particular problem does not occur in 
field theory, there remain severe difficulties in carrying out the programme. 
The most promising approach (3) put forward in Ref. [S| still contains un- 
cancelled Jacobian factors, which can only be represented in the Lagrangian 
by the introduction of several auxiliary fields. 
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